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The growth of supermassive black holes in the centers of galaxies: 
absorption of stars and activity of galactic nuclei. 

Abstract 

The evolution of a supermassive black hole in the center of galaxy is consid- 
ered. We analyze the kinetic equation describing relaxation processes associated 
with stellar encounters. The initial distribution function of stars is assumed to be 
isothermal, and it changes in accordance with consumption of stars by the black 
hole. We calculate the consumption rate and obtain the law of growth of super- 
massive black hole, i.e. the dependence of the black hole mass on time and on 
the velocity dispersion of stars. The obtained values of mean black hole mass and 
mean rate of stars consumption by a black hole are in reasonable agreement with 
observational data. 

1. Introduction. 

The massive compact concentrations of mass in the centers of galaxies are believed 
to be supermassive black holes with masses 10 6 -j- 10 9 solar mass (M ). The recent 
q ■ observations [I] provide significant evidences for this hypothesis. 

There are two global approaches to the origin of these black holes: they could either 
form together with the galaxy or even earlier j2j, or grow up during galaxy evolution, 
capturing matter from the surrounding. We throughout the paper consider the second 
approach. 

Up to the present time there has been no consistent theory how the growth of black 
holes in the galactic centers goes on. 

Black hole can capture luminous (stars and gas) and cold dark matter from the central 
parts of galactic core (bulge). Consumption rate of stars has been investigated in jHlSE! 
EMU. (The first statement of the problem of distribution and absorption of particles in the 
Coulomb potential has been considered in [Hj). In [3101111 the diffusion in two-dimensional 
energy-momentum space has first been discussed as the main process. However, all these 
authors considered the dynamics of stars in the Coulomb potential of black hole, whereas, 
as it has already been pointed out in jHl , in the case of giant black holes the consumption 
rate is determined by distances where it is negligible and the stars move in their own self- 
consistent potential field. On the other hand, in all these papers stationary solution has 
been considered; as we will show below, the solution of the real problem is non-stationary. 

The kinetic theory describing dynamics of non-baryon cold dark matter (CDM) alone 
(without baryon fraction and central black hole) has been developed in |9 3 . It has been 
shown that the dark matter should have formed quasi-stationary spherical structures with 
singular density distribution, bound by gravitational forces. These structures are now 
observed as galactic halos. The baryon matter should have fallen after recombination in 
the potential wells formed by dark matter, beginning formation of galaxies and forming 
seed black holes in the centers of the wells jlOj . This picture explains well why the black 
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holes are situated exactly in the dynamical centers of galaxies the potential is very 
smooth in the central part, and the black hole could not reach the center unless it would 
be formed just in the center, in the potential well already existing. 

The contemporary infall of CDM particles into the black hole caused by gravitational 
interaction with stars and the corresponding growth of black hole has been investigated 
in • It was shown that this process leads to formation of black holes whose masses are 
consistent with the most part of supermassive black hole masses measured nowadays. 

In the present paper we analyze the mutual dynamics of cold dark matter and stars 
in central parts of a galaxy. In accordance with previous papers, we find that the most 
important process leading to the black hole growth is (for both components) far encounters 
causing the relaxation in the energy-momentum space, or, to be more precise, in the space 
of adiabatic invariants. The influence of black hole appears both in gravitational potential 
(i.e. the view of adiabatic invariants) and - even stronger - in boundary condition. 

We derive the equations describing time- dependence of stars- and dark matter- distri- 
bution functions, and calculate the coefficients in one particular assumption of isothermal 
distribution of stars and the potential dominated by stars, either (Section 2). 

In Section 3 we discuss the boundary condition caused by presence of black hole, and 
the influence of the black hole potential. 

Further, we find an approximate solution of the equation for stars, ignoring presence 
of dark matter and gas (Section 4). The result shows that black hole mass should grow in 
time as t 1//2 , and is proportional to a 3 / 2 , where a is the rms velocity in the inner regions 
of bulge. 

In discussion, we review briefly the results of the paper. We show that the mean energy 
released due to stellar matter absorption is in reasonable agreement with the observed 
AGN activity. 

2. Two-component dynamics. 

We denote by /*(r, v) and /<i(r, v) the distribution functions of stars and dark matter, 
respectively. This means that fraction of mass of sort a contained inside the interval d 3 r, 
d 3 v is f a d 3 rd 3 v. 

The mutual dynamics of dark matter particles and stars including diffusion in the 
momentum space is described by kinetic equation with Landau scattering integral in the 
right part [T3]: 



or 




Here greek indexes denote the sort of particles, M a denotes the mass of particles, f a (r, v) 
corresponds to scattered particles, /^(r, v') corresponds to scatterers, 

w kp = (u 2 5 kp - u k u^j /u 3 , u = v' - v , 
u is the relative velocity, A ~ 15 is gravitational Coulomb logarithm. The gravitational 



2 



potential \I/(r) is determined by Poisson law, 

A^ = 47rG5> Q (r)=47rG]r J /„(r,v) d 



3 v 



where p Q (r) is mass density of fraction a. 

Comparing the terms in ((TJ one can see that scattering of both dark matter and stars 
by stars is M*/Md m times more effective than scattering by dark matter particles, and in 
the case of dark matter the second term is M<2 TO /M* times smaller than the first. So we 
rewrite the equation JI} in the form 

^ = 2vrG 2 M«A (^-W^ - £-W k f.) (2) 
dt \ov k ovp ov k J 

^ = 2nG 2 M^W kp ^ (3) 
at ov k ovp 

where 

W k = J ^w kp d 3 v' ,W kp = J f*w kp d 3 v' , (4) 
Rewriting equation in terms of adiabatic invariants 

If the potential is spherically symmetric, one can introduce Hamiltonian variables 
(1,0), where I = (I,m,m z ) are adiabatic invariants (m is absolute value of angular 
momentum, m z is its projection on z axis, I is radial action; all - per unit of mass): 

r+(E,m) 

v r dr 



2.7--, J f (5) 

r_(E,m) ' 7 
2 2 

E = — + *(r) , m = r x v , £? - *(r±) - —f = , 

and 4> are corresponding angles. One of the advantages of these variables is that the 
adiabatic invariant / of a particle (unlike the energy E) does not change during slow 
changes of the potential. 

In these variables, the left side of the kinetic equation does not contain derivatives 
over action variables I k [T^ ITi] : 

d/q(I» m df a _ , , 

We now make use of the fact that the frequency of collisions is much smaller than the 
frequency of orbital motion. After averaging over the angle variables <p k the second term 
in the equation vanishes, and the equations then take the form (see [TH H2] for 

detailed derivation) 

dfJX) in dU_M^^_ 
dt dl k kp dl p M*dI k KkTa} ' [) 
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where 

1 r dT 



P ' 

(7) 



j3 dl k dl p 
(27r) 3 J dv m dv ri 



R kp = 27rG 2 M Ht A T ^^ / d 3 (t> „ „ P W„ 



We assume that in the vicinity of the center the distribution function is spherically sym- 
metric, i.e. does not depend on the direction of angular momentum. Then, there is a 
cylindrical symmetry in the space of adiabatic invariants. Making use of 

-K3 — n 2 > -n-31 — -n-21 > -"-32 — -"-22 , 

m m m 

we can exclude m z from the equation (jUJ). Eventually, we obtain the equations for distri- 
bution functions of stars and of dark matter: 

at ol \ a I am J m am \ al am J 

— = — (r — + R —\+——m(R — + R —\ (9) 

dt dl \ 11 dl 12 dm J m dm \ 21 dl 22 dm J 

We see that evolution of distribution functions is described by axially symmetric diffusion 
in three-dimensional space of adiabatic invariants. 

We note that equations (JHJ), © have a divergent form, i.e. 

- = —s k 

dt dlk 

The flux of matter through any surface in the momenta space is given by the integral of 
S over the surface. For example, the flux through the plane / = is J dm J dm z Si = 
J 2mSidm. The boundary / = is not physical, hence flux of matter over any part of it 
should be zero. Thus, we obtain the boundary condition 



Z>i\i=o - [ Kll ~ J~ + 12 5^ ~ Kl * a 



= (10) 

7=0 



In absence of a black hole, the other boundary condition would be 

The second boundary condition imposed by presence of black hole will be considered in 
Section 3. 

Isothermal distribution function 

One important example is the Maxwellian stars distribution function with arbitrary 
time-independent potential: 

f M = foe- E/ ° 2 
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This function corresponds to the equilibrium state, since it satisfies equation (J2J). This 
means that both parts of equation (jSJ) are equal to zero. Moreover, flux of matter over 
any surface in the momenta space must be zero, hence all components of the flux S must 
vanish. This puts two relations on coefficients Rk,Rk P - 

7? ^ i p pf _n z? ^ i p of _n 

-n-n-^r- + iii2T; ttijM — U , ^121-?:^" + ^22 t; -K2/M — U (^iij 

oi am a7 am 

In the particular case of potential dominated by stars, there is the only parameter a 
characterizing the distribution function: 

fis = foe~ E/a2 , * is = 2a 2 In r , f = [(27r) 5 / 2 Ga] ~* . (12) 

This function is called isothermal, the parameter a 2 corresponds to velocity dispersion. 

We now calculate the coefficients Rk, Rk P for isothermal distribution function, /* = f is . 
All manipulations needed for the procedure are described in Appendix. The result is 

R 22 ~ 0.5GM*A(T = R 
R u ~ -0.6i2(l - //$ a (/i)) 
R u ~ fl(0.4$i(^) - 0.8/i$ 2 (/i) + 0.4) 

fli - - / + Q 6m (0-7$i(/i) - 0.8^2^)) ; (13) 

R > * - L3 7T^^ 2(/i) ; 

m 

/i ~ 0.3- 



/ + 0.6m 

where approximate expressions for functions $ are: 

$ 2 (//) ~ 1 + 0.35fi~ 06 , 

$x(/i) ~ 0.5 - 0.8 ln/i • ($ 2 (/i) - 1) ~ 0.5 + 0.33^ a9 . 

The isothermal distribution function can also be expressed in terms of Hamiltonian 
variables (see Appendix): 

/,(/,-») =/„e-"^^g, (14) 



3. The black hole, and setting the problem. 

The scattering of dark matter particles by stars and their consumption by a black 
hole has been discussed in |12j . It has been assumed in that paper that growth of black 
holes is fully determined by dark matter. Here we discuss an opposite assumption of full 
absence of dark matter. From now on, we discuss only dynamics of stars. Accordingly, 
we hereafter omit the index * in 
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Black hole can effect the dynamical evolution of surrounding matter in two ways: 
consuming matter from the loss-cone and changing the gravitational potential. We now 
consider the first effect as principal one, and later on will briefly discuss the second. 

The black hole with mass Mbh disrupts (and after that consumes) all stars that come 

nearer than the tidal radius r t = (f^^ 1 ) , where is the density of the star. We 
assume it to be approximately equal to solar density p . In terms of variables (I, m) this 
means that during one oscillation period all stars inside the loss-cone 

/ hjt \ 2 /3 / \ 1/6 
as / M bh \ Pe \ _2 



m < m t = pGM bh r t ~ 6 • 10 18 1^1 cm 2 /s (15) 

are disrupted. This boundary condition on m has first been pointed out in j^J. On the 
other hand, if angular momentum is less than m g : 

m<m g = AGM bh /c , (16) 

the star will anyway fall into the black hole [T^], independently of m*. (If the black hole 
mass Mbh is greater than ~ 3 ■ 1O 8 M , the tidal radius becomes smaller than the gravita- 
tional radius r g , and stars become consumed without disruption. For dark matter particles 
the loss cone is always defined by m g .) For brevity, we shall hereafter denote the boundary 
value by m t , assuming maximal of the two values m g and m t : m t = max{m 9 ,m t }. 

One can easily estimate the growth of the black hole caused by direct capture of matter 
from the loss-cone: taking, e.g., the isothermal function ()14[) . we have AMbh ~ 2iramt/G. 
For black holes with masses more than 1O 3 M relative increase of mass AMbh/Mbh is 
less than 10% . Hence the main growth is connected with diffusion of matter into the 
loss-cone. 

The condition that the particles vanish as they get into the area m < mt may be 
expressed in terms of equation (jHJ) by boundary condition 

f\m=m t = (17) 

It is the change of gradient of distribution function df /dm near the boundary that causes 
shift of particles in the momenta space into the loss-cone. This mechanism may lead to 
consuming much more mass than has initially been inside the loss-cone. 

Applying Gauss theorem to the equations ©,© we see that the rate of black hole 
growth dMbh/dt, i.e. loss of whole mass of stars in the surrounding space, is equal to flux 
through the boundary m t : 



+m t 



-ra t 



M bh = S= J d 3 jdl J dm z S m = {2itf2m t J dl (#22^ + i?i 2 |j - R 2 f 







m=mt 

fl8) 



The last two terms are zero because of the boundary condition (117(1 . 

So, to find the contemporary mass of black hole and distribution of stars in its vicinity 
one should solve equation (jHj) with boundary conditions (jlOJl and ()17j) . The isothermal 
function (|14|) is the most natural choice of initial distribution because it is stationary 
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distribution with self-consistent potential and with no flows. If there were no black hole, 
the isothermal distribution remained constant. The boundary condition (|17|) disturbs the 
equilibrium and makes the distribution change. 

Furthermore, we take the same isothermal function to calculate the coefficients Rk, Rkp 
at any moment of time. We suppose that one may neglect the changes in distribution 
of scatterers, because (J7|) includes only integrated function and changes in do not 
effect the coefficients significantly. We will check the convenience of this approximation 
in Section 4. 

Let us now estimate the influence of gravitational potential of black hole. Taking 
values Rik of the isotherm ()13|) for an estimate, one can see that to the lifetime of the 
Universe t ~ 3 • 10 17 s different members of (jHJ) may produce significant changes of / at 
scales m ~ / ~ y/Rt ~ 8 ■ 10 25 cm 2 /s , this corresponds to 



Dynamics of matter at these distances determines contemporary infall into the black hole. 
This is much smaller than size of bulge (~ lkpc), so we should not take boundaries of the 
distribution into account. On the other hand, the black hole potential becomes significant 
at distances 



the region of influence of the black hole potential does not reach up to characteristic 
distances (119)) . In this case only a small fraction of particles enters the region, and, 
because of ellipticity of most orbits, spends there rather small part of time. This is in full 
conformity with conclusions by j^j. 

The same arguments are applicable to the region r < r co u ~ R Q Mbh/M Q where orbital 
velocities become comparable to the escape velocity from a typical star, and inelastic 
scatterings are important E] • One can easily check that r coU ~ r a for values a under 
discussion. 

4. Asymptotic solution. 

The equation (JHJ) is difficult to solve. In this section we obtain at least some integrated 
characteristics that make us able to estimate the total consumption rate S of stars and 
its evolution in time. 

Let us first discuss the vicinity of the black hole boundary, m t < m <C \fR^t and 
I > in. Comparing different terms in equation (jHJ) we find out that, e.g., the term 
containing i? 2 2 is of the order of -R22/A™ 2 ^ df /dt, which means that the picture near 
the boundary is quasi- stationary at contemporary time. All the other terms except the 
two (with Ru,Ri) that are singular at m — > (see are maximum of the order of 




(19) 




We see that for not very massive black holes, 




(20) 
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Risf /I™, i- e - much smaller. To get rid of the two remaining terms, we integrate (JBJ) over 
I and introduce F(m) = J f(l,m)dl. Thus, we obtain 1 

mR 22 ^— ~ (21) 



m dm dm 

Taking into account the boundary condition at m = m t , we find 

F(m) = C(t) In — , m < VRt (22) 
m t 

The unperturbed isothermal function (fT4*|) gives a solution far from the black hole bound- 
ary, 

F(m) = 0.5f a 2 — , m > y/Rt 
m 

Substituting (|2*2*Jl into (fTHjl. we see that 

5 = (27r) 3 2 J R 22 C(t) (23) 
The constraint on C(i) is that the flux S is equal to the whole loss of mass: 

S = J d 3 <p J dl J dm J dm z f = -(2n) 3 —J 2mF{m)dm 

—m 

Comparing this with (J23j) and calculating the integral, assuming the two asymptotes are 
equal at m — mi, we obtain two equations: 

C(t) mi \n— = 0.5f a 2 , 

m t 

C V m t 2/ 

If the boundary does not move, m t = const, then the approximate solution to the problem 
is 



(t) ~ 2v / ^ln" 1 / 2 f^E] 



m x {t) ~ 2V J Rtln" i/2 — (24) 



C(t)^^Lln- 1 / 2 f^') (25) 

Now, knowing the dependence of flux S = M^h on time, we calculate the expression for 
mass of black hole. 



Lbh 



(2nYf aWRt 



\ 



, 2VRt 
m 

m t 



1 We neglect here small values of J, for which the terms with R12, R2 may also be significant (although, 
no more than of the same order as the term with i?22)- However, it is easy to check that in the vicinity 
of the boundary (and hence, everywhere) the adiabatic invariant of each star could only grow with time 
(i.e., Si < 0), and thus, contribution of small values / in the flux S and in the integrated function F is 
not large. 



S 



Substituting value of / , we find 



M bh ~ f J= (26) 



mi 



The remarkable fact is that the result depends only weakly on the boundary m t . This is 
why the solution (|24 |) -(|26 |) remains the same (up to the same accuracy) if we account the 
change of m t with time ((jl5j) or ([16))) according to change of the black hole mass (|26|) . 

Taking A = 15, 3 • 10 17 s - the Universe lifetime, and mass of star M* equal to solar 
mass, we obtain the mass of black hole at arbitrary time: 



3/2 



M^1.2.1O'M fex5^j i¥10^ (27) 
The contemporary consumption rate of stars by black hole is (|23j) 

S « ^ » 6 • l(T 4 M /year • f ^— ^1 ^ . (28) 

2t 0/ V200km/ S ; v 7 

It decreases in time as t -1 / 2 , and corresponds to 3.6 • 10 43 erg/s at contemporary time if 
a = 200km/s. Occasionally, at a ~ 200km/s tidal and gravitational limits of momentum 
mt and m g corresponding to the mass (|27J) are close to each other, 

m g ~ 2.7 ■ 10~ 3 VRt nnr ? , , m t ~ 4 

The corresponding mi is 



v '""^OOkm/s ' *" v ""V 200km/s 




mi ~ 0.8Vi?t . 

We note that mass of black hole is no more than half the mass that was initially (when 
there was isothermal distribution) in the area m < mi. The same is correct for any 
time (excluding the very beginning). This proves our initial assumption that change of 
distribution function does not affect significantly the coefficients Rk p - 



Possible constraint. 

Let us now discuss briefly possible constraints on diffusion model. All equations begin- 
ning with (jUJ) are averaged over orbital period. The boundary condition (|T7|) also assumes 
that all particles with average angular momentum less than mt become captured into the 
black hole. This means that during one period the change of trajectory is negligible. But 
for very elongated and large trajectories change of angular momentum of particle dur- 
ing one orbital period is comparable with the boundary angular momentum m t . Hence, 
whether the star will or not be captured by black hole during one period depends not on 
the mean trajectory, but on its small fluctuations caused by individual collisions. 2 In 
terms of adiabatic invariants this may be formulated as follows. 

2 The effect was first pointed out in 0J|S]. It was then formulated in terms of critical energy, or critical 
radius in real space. 
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The diffusion approximation and equation (JSJ) are applicable at any point of the mo- 
menta space. However, the boundary condition (|17|) may not be applied at very large 
values I (i.e., for very long orbits). Namely, mean change of angular momentum Am 
during one period T(J, m) is of the order of Am ~ yRJ^J^m), hence the constraint on 
CZJ is 

yjKT(I,m t ) < \J RT(I crit) m t ) = m t 

On the other hand, we have seen in Section 3 that main changes in the distribution 
function, and hence main contribution to the flux occur at scales I < \f~Rt where t is the 
age of the black hole. This means that the consumption rate of stars is not disturbed by 
influence of l cr a if 

/ ~ \f~Rt < I crit 

At sufficiently large values I ^> m t the orbital period is approximately equal to 
T(I,m t ) ~ i\l jo 2 (see Appendix). Substituting m t from (fTo} we obtain that one anyway 
should not take the existence of I cr u into account if the black hole mass is greater than 

/ 1 t 200km/s 

M bh > - — 

- \50 3- 10 17 s a 

Nowadays, this corresponds to minimal ~ 2- 10 6 M Q . Comparing the time dependence 
in ()29|) with that in (|2Tj) we find that (if o has not changed dramatically) the masses of 
the black holes that have grown as a result of the described relaxation mechanism have 
always been greater than the critical limit as far as galaxies exist. On the other hand, for 
black holes with sufficiently smaller masses additional investigations of influence produced 
by Icrit are needed. 







(29) 



5. Discussion 

In the paper we derived the equations describing kinetic relaxation of two-component 
(stars + dark matter) system in the vicinity of supermassive black hole, in self-consistent 
gravitational potential. The equations have been expressed in terms of Hamiltonian vari- 
ables. 

We have analyzed the growth of black hole caused by consumption of stars, without 
taking non-baryon dark matter and luminous gas into account. We have seen that for 
black holes smaller than 1O 6 M (J29j) individual collisions near pericentre might play crucial 
role. For black holes larger than 10 8 M Q ()20|) one should take the potential of the black 
hole into account. However, the most part of discovered black holes JT] have masses in 
the range between 1O 6 M and 1O 8 M where the main mechanism leading to consumption 
of stars is diffusion in the momenta space, or in the space of adiabatic invariants. 

Analogous mechanism has been discussed previously in a set of papers [13 IH El El E| • 
But all of them have examined stationary distribution in the Coulomb potential of the 
black hole. We have considered the evolution in self-consistent potential and found that 
there is no stationary state in presence of black hole. Also in previous papers only diffusion 
along one axis (either m or E) was taken into account. The advantage of Hamiltonian 
variables m, / is that in this variables the diffusion along m axis contains both shifts 
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along m- and EFaxes in the (E, m) plane and makes the main contribution into the flux 
on the black hole. The diffusion along I axis plays less important role in the absorption 
process, being merely responsible for change of form of stars trajectories. Flux along I- 
axis increases smoothly as m decreases, reaching its maximum at the black hole boundary 
m = m t , where Si/S m = —0.6 (i.e., the particles in the I — m plane intersect the black 
hole boundary m = m t at the same angle independently on /). 

The solution (|2*7|) obtained in Section 4 gives reasonable estimate for contemporary 
characteristic black hole mass which lies in the middle of observed range of masses. The 
estimate for flux of matter onto a black hole (J2"H)) also corresponds to average observed 
activity of galactic centers. 

Equation (J2*7)l also shows that M bh should change very slowly in time, M bh ~ \ft ~ 
(1 + z)~ 3 ^. Hence, the corresponding difference between black hole masses in far and 
nearby galaxies is too small and hidden by other effects. On the other hand, the obtained 
M — a relation differs from that derived from observations (from ~ a 3 7 in jTHj to ~ a 5 3 
in [TZj) However, this discrepancy is not critical, since black hole masses depend not on 
one parameter a only, but also on the bulge mass ^1] , mass of halo [T%1 IT2*] , amount 
of interstellar gas etc. For example, the contribution of dark matter could increase the 
resulting black hole mass in those galaxies where amount of dark matter in the central 
parts is significant. 

In the previous paper [TJ] we have discussed the growth of black hole caused by 
consumption of dark matter scattered by stars. The result was that the characteristic 
black hole mass is 

19 ^ 4 

M mM) = 8 • IV M (^) ' (J""-) " ' (^) 1 , (30) 

where Mjj and Rh are the mass and the radius of the dark matter halo. 

This result and solution (|2*7j) correspond to two limit cases describing the black hole 
growth caused either by dark or by luminous matter consumption. We see that contri- 
butions of stars and dark matter into the black hole growth are comparable for galaxies 
where densities of dark and luminous matter in the bulge are of the same order. But dark 
matter plays in this case a dominant role. This, in particular, means that in galaxies 
where observations show very small amount of dark matter the black hole mass should be 
in average smaller than in others. With account of real variations of all parameters (J27|) 
and (|3L)p cover the range 

5 • 1O 6 M <M bh <2- 1O 8 M (31) 

which includes the most part of observed black holes. 

The common point of view is that the consumption of stars by a black hole in the 
galactic centre is the main source of energy providing activity of the galactic nuclei [T^ I2Uj. 
To verify this hypothesis we compared the observed X-ray luminosity of the galactic 
centres with the central black hole masses for 46 galaxies (Table 1). From the Table it 
follows that: 

1. The most part of the black holes have masses in the range (j3~Ij) restricted by ()27|) 
and flU. 
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2. The observed intensity of X-ray emission is in a reasonable agreement with the mean 
energy released as a result of stars consumption (J25|) . 

3. The analysis of the data shows very strong variations in observed luminosity of galac- 
tic nuclei. But no correlation with the black hole masses was found. This indicates 
that the luminosity of one nucleus should vary strongly in time. In other words, 
the duration r e of energy emission following a star absorption is much less than the 
mean time interval between two absorption events T ~ M Q /S ~ 2000 years (|2~%j) . 
From observed variations (3-4 orders of magnitude), assuming an exponential fall 
of emission in time, we find an estimate r e < T/10 ~ 200years. It follows that the 
most part of galaxies containing supermassive black holes are not active at present 
(as, e.g., Milky Way). On the other hand, maximal intensity of emission caused by 
star consumption should be much more than the average value 5* obtained in ()28|) . 

4. In Fig.l, the histogram of average luminosity versus the black hole mass is presented. 
The histogram is based on Table 1. It shows clearly that the luminosity breaks 
sharply at M&/, > 3 • 10 8 M Q . This corresponds to the fact mentioned in Section 3 
that for such black holes the gravitational radius is greater than the tidal radius, 
and stars sink into the black hole without disruption and energy emission. The 
observed cutoff proves that it is the infall of stars - not gas etc. - that provides the 
main contribution in the emission activity. On the contrary, in quasars with large 
redshifts the main mechanism causing energy emission may be consumption of gas 

m 

Thus, in this paper and in we present the theory of growth of supermassive black 
holes in the galactic centers caused by Coulomb scattering on stars and further consump- 
tion of stars and dark matter. The theory is in reasonable agreement with observations. 
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Table 1. 

The X-ray luminosities and the central black hole masses of galactic nuclei. 3 



Name 


L, ere 


;/s 


M bh , M a 

UIL7 


3C 120 (Mrk 1506) 


9.772 


•10 43 


2.3 • 10 7 




1.343 


• 10 44 




Ark 120 (Mrk 1095) 


7.805 


•10 43 


1.84- 10 8 


Circinus 


2.47- 


10 40 


1.3 ■ 10 6 


Fairall 9 


6.688 


•10 43 


8.0 ■ 10 7 


IC 1459 


8.398 


.1 40 


3.7- 10 8 




1.051 


•10 41 




IC 4329A 


3.748 


•10 43 


5- 10 6 




2.541 


•10 43 




UGC 3973 (Mrk 79) 


2.594 


•10 43 


5.2- 10 7 


Mrk 110 


1.625 


• 10 44 


5.6 • 10 6 


Mrk 335 


2.04- 


10 43 


6.3 • 10 6 




4.916 


•10 43 




Mrk 509 


1.188 


• 10 44 


5.78 • 10 7 


Mrk 590 (NGC 863) 


3.853 


•10 43 


1.78 • 10 7 




8.884 


•10 43 




NGC 205 (MHO) 


< 1.172 • 10 38 


< 9.3 ■ 10 


NGC 598 M33) 


1.46 • 


10 39 


< 1.5- 10 




1.902 


.1 39 




NGC 1068 (M77) 


1.315 


•10 41 


1.6 • 10 7 




5.124 


•10 41 






8.81 • 


10 41 




NGC 3115 


1.773 


.1 39 


9.1 • 10 8 


NGC 3227 


1.371 


•10 42 


3.9- 10 7 




7.127 


•10 41 




NGC 3516 


5.578 


•10 41 


2.3 ■ 10 7 




1.084 


.1 43 




NGC 3608 


1.174 


.1 40 


1.1 ■ 10 8 


NGC 3783 


8.515 


.1 42 


9.4 • 10 6 




7.384 


•10 42 




NGC 3998 


1.425 


•10 41 


5.6 • 10 8 


NGC 4051 


7.609 


•10 41 


1.3 • 10 6 




7.624 


•10 41 




NGC 4151 


9.025 


•10 42 


1.53 • 10 7 




5.578 


•10 41 






4.784 


•10 41 





3 The black hole masses are adapted from |TT], the X-ray luminosities - from |21|. Only galaxies with 
measured X-ray luminosities are taken from the sample of 
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NGC 4203 


7 


471 


• 10 40 


< 1. 


2 • 10 




6 


424 


• 10 40 






NGC 4258 (M106) 


2 


914 


• 10 40 


4.1 • 


10 7 


NGC 4261 (3C 270) 


1 


575 


• 10 41 


5.2 • 


10 8 




1 


434 


• 10 41 






NGC 4291 


6 


128 


i A(\ 

■ 10 40 


1.5 • 


10 s 




6 


769 


• 10 40 








4 


28 • 


10 40 






NGC 4342 


2 


077 


• 10 39 


3.4 • 


10 8 


NGC 4374 (M84) 


6 


317 


• 10 40 


1.6 • 


10 9 




6 


657 


• 10 40 








6 


138 


■ 10 40 






NGC 4459 


4 


169 


• 10 39 


6.5 • 


10 7 


NGC 4473 


1 


109 


■ 10 4U 








5 


089 


• 10 39 


1.0 • 


10 8 


NGC 4486 (M87) 


3 


257 


• 10 42 


3.4 • 


10 9 


NGC 4593 


6 


535 


• 10 42 


8.1 • 


10 6 


NGC 4594 (M104) 


3 


361 


• 10 40 


1.1 • 


10 9 


NGC 4649 


1 


002 


• 10 41 


2.0 • 


10 9 




1 


587 


• 10 41 






NGC 4697 


9 


141 


• 10 39 


1.2 • 


10 8 


NGC 4945 


5 


447 


• 10 39 


1.1 • 


10 6 


NGC 5548 


2 


182 


• 10 43 


1.23 


• 10 s 




2 


154 


• 10 43 








2 


778 


• 10 43 






NGC 6251 


1 


71 • 


10 42 


5.4 • 


10 8 


NGC 7469 


1 


699 


• 10 43 


6.5 • 


10 6 




3 


34 • 


10 43 








2 


071 


• 10 43 






PG 0026+129 


2 


798 


• 10 44 


4.5 • 




PG 0052+251 


4 


766 


• 10 44 


2.2 • 


10 s 


PG 1211+143 


1 


098 


■ 10 


4.05 


1 n7 

• 10' 


PG 1411+442 


5 


049 


•10 42 


8.0- 


10 7 


PG 1426+015(Mrk 1383) 


8 


507 


•10 43 


4.7- 


10 8 


PG 1613+658(Mrk 876) 


1 


652 


• 10 44 


2.41 


• 10 8 


PG 1617+175(Mrk 877) 


6 


188 


•10 43 


2.73 


• 10 s 



Appendix. 

Calculating the coefficients R k ,R kp in ((7j) and the adiabatic invariant in the 
case of isothermal stellar distribution function. 

A number of simplifications may be made for any isotropic distribution function. We 
first discuss the coefficients Wk p ,Wk introduced in (j3J). If the distribution function /* is 
isotropic at any point of coordinate space, i.e. does not depend on momentum, but only 
on energy, then tensor W kp depends on only one vector v, and hence takes the form 

W kp = A(E,r)5 kp - B(E,r)^, (32) 

where, calculating convolution W kk and product W kp v p , we find 

°? ( % „> I ,,'2 



A=* [ dE >f(E>)lii(l-m . E'< E 
3 J JK M 1 , E' > E 

oo 

A - B = ^ / dE>f(E> 



oo r n 

1 , E' > E 



j2(E-*(r)) , v' = yj2(E>-*(r)) 



and similarly 



W fc = ^D(£,r), (33) 



oo 

87T f 7 771' 9 



To calculate the coefficients R k ,R kp in (jZJ it is useful to introduce new variables £ 
{CiiG} = {E,m}. Then, since Jj = Ii(E,m), we have from (JJJ) 

For brevity, we denote 

An = 2ttG 2 M*A , < .. >= — — / 



(2tt) 3 

Then, making use of (J32|) and we obtain the expressions for R' a , R' ab : 

R[ = An < vD > ; i?2 = < D/f > ; 

(34) 

R' u = An < w 2 (A - 5) > ; R' 12 = A m < A - B > ; R' 22 = A < Ar 2 - Bm 2 /v 2 > 

Note that because of spherical symmetry of all coefficients, only one integral remains in 
the averaging: 



2\E- *!ri-?4 



T(i?, m) J v r ' ' / f r ' \ \ 2r 
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All the previous results are applicable for any isotropic function. Now we substi- 
tute the isothermal distribution function (|12j) . For calculations it is useful to introduce 
dimensionless variables: 



"' x = re- E ^ 2 



2oe E l 2 ° 

The condition v 2 = 2(E - * - m 2 /2r 2 ) > leads to 

< x < 1 , 0</i< \j\f2e . 
Then, for example, period of oscillations is 

T(E, m) = V2[ d ^E-*-™ 2 /^) = l^f^ ? f{fi) = j dxOj-lnx - , 2 / x 2 ) 
J yJE - m? /2r 2 ° { ^-lnx-fi 2 /x 2 

The remarkable fact is that T(fi) is almost constant! It grows very slowly from 1.77 at 
jj, = to 1.9 at /i ~ 0.4 - maximum possible \x. 

The averaging procedure now can be rewritten as 

2 r dr 1 r dx9(— In x — [i 2 /x 2 ) 

T{E,m)J_ V r [ - ] ~TiiT)J o ^-\nx-fi 2 /x 2 

We also can calculate A, B and D for the isothermal distribution function: 

A=— f -±o 2 (-2 \nx(2H + x 2 ) - H + x 2 ) , A - B = — fya 2 (2H + x : 



1 



Z^ = _|^(2tf + x 2 ) , W 2 = -4a 2 lnx, H{x) = J ydy^ 

X 

Note that 

A - B = -a 2 D/v 

which leads to (jllj) . 

From (|34|) we now find 

R'xx = — ct 2 R'i , R'12 = — R*2 
R'i(E, m) = -4a 2 e-^ 2 A y/o$i(^) 
R' 2 (E,m) = - me -^ 2 A y/o$2(/i) 

87T 

R' 22 {E,m) = a 2 A y/o$2 2 (/u) 

where 

^> 2 (aa) =< 2H(x)/x 2 > +1 , =< -Inx ■ 2H(x)/x 2 > + < - \nx > 
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$ 22 ( /i ) =< -2\nx{2H + x 2 ) - H + x 2 > -fx 2 < AH/x 2 + 1 + 3H/(x 2 lux) > 
Up to at least < 10~ 4 one can use approximate expressions: 

& 2 (jJ,) ~ 1 + 0.35/i- a6 , 

~0.5-0.81n/i-($ 2 (/i)-l) ~ 0.5 -0.3 In // ■ // _a6 , 
$ 22 (/i) ~ 0.87 
Finally, we calculate the adiabatic invariant I 

/(£, m) = ^ e E/2a ' 2 J ^/-\nx- fi 2 /x 2 6(- In x - fi 2 /x 2 )dx 



o 



Calculating the integral we obtain 



/(£, m ) ~ ^ e ^ 2 (o.9 - 2/x) = -(0.9ae E ^ 2 - m) 

IX IX 



Hence 



2 „ 7 + -m 



m ~ ix ' 2a 2 ' 

The distribution function can now be expressed in terms of adiabatic invariants: 

'2-0.9\ 2 1 



tx ) (I H — m 

v 7T 

Another useful formula is the expression for the period T(E, m) formulated in terms of 
J, m: 



T(I,m) 



a a 0.9a 
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